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We study transport in two-terminal metal/quantum spin-Hall insulator (QSHI)/metal junctions. 
We show that the conductance signals originating from the bulk and the edge contributions are not 
additive. While for a long junction the transport is determined by the edge states contribution, for 
a short junction, the conductance signal is built from both bulk and edge states in the ratio which 
depends on the width of the sample. Further, in the topological insulator regime the conductance for 
short junctions shows a non-monotonic behavior as a function of the sample length. Surprisingly this 
non-monotonic behavior of conductance can be traced to the formation of an effectively propagating 
solution which is robust against scalar disorder. Our predictions should be experimentally verifiable 
in HgTe QWs and Bi2Se3 thin films. 



I. INTRODUCTION 

Topological states of mat ter a re characterized by bulk 
invariants: Chern numbers^ or by a Z2 invariant 3 , 
depending if the time reversal symmetry is broken or 
conserved. Due to the bulk-boundary correspondence^^ 
this classification can be translated to the existence of 
topologically protected states at the edges of the system. 
In the case of the Z 2 insulator the edge states are formed 
by time reversed modes, so called Kramers' partners 
and they are helicaP, i.e. Kramers' partners counter- 
propagate along a given edge of the sample. Depending 
if the number of Kramers' pairs on the edge is even or 
odd, the system is topologically trivial or non-trivial. 
Although the topological invariant is the property of the 
bulk, experiments usually detect edge states properties. 
Indeed the confirmation that HgTe/CdTe quantum wells 
(QWs) are two-dimensional topological insulators^ was 
provided through the measurement of quantized edge 
conductance^"^. 

One benefit of HgTe/CdTe QWs is that the topologi- 
cal order can be controlled by the thickness of the HgTe 
layer: below the critical thickness of d c ps 6.3 nm the sys- 
tem is a trivial insulator, whereas above d c the system 
behaves as topological insulator. This gives the possibil- 
ity to observe the topological phase transition directly. 
Therefore it is of great interest to find further experi- 
mental measurable indicators which distinguish the two 
regimes. Since the topological invariant is a bulk prop- 
erty, it seems natural that the bulk conductivity could 
also carry information about topological properties of 
the system. Indeed it was proposed in Ref. [T2] that a 
non-monotonic conductance as a function of geometri- 
cal aspect ratio in metal/HgTe/metal junctions with a 
characteristic maximum describes the topological insula- 
tor regime. The conductance maximum occurs when in- 
coming metallic bulk states tunnel through a short and 
wide two-dimensional topological insulator and is robust 
against scalar disorder on the order of the energy gap^. 
However, so far this analysis was limited to the periodic 



boundary conditions (PBC) and neglected the existence 
of edge states, which can significantly modify the behav- 
ior in experimentally relevant setups. Here we analyze 
carefully the properties of the single QSHI/mctal inter- 
face as well as the interplay between edge states and bulk 
states in metal/QSHI/metal junctions. To do so we ap- 
ply a generalized wave matching method based on Fourier 
modes, like it was e.g. used to analyze two interface tun- 
neling junctions in HgTe 14 . Zhang et al. restricted them- 
selves to a low energetic behavior around the gap, which 
is dominated by the linear dispersion of the bulk and the 
edge states. In contrast our studies combine a topologi- 
cal insulator with two highly doped metallic leads treated 
as highly doped HgTe with quadratic dispersion. We 
show that the conductance signals originating from the 
bulk and the edge contributions are not additive. While 
for a long junction the transport is determined by the 
edge states contribution, for a short junction, the con- 
ductance signal is built from both bulk and edge states 
in the ratio which depends on the width of the sample. 
Further, the conductance for short junctions shows a non- 
monotonic behavior as a function of the sample length in 
the topological insulator regime. Surprisingly this non- 
monotonic behavior of conductance can be traced to the 
formation of an effectively propagating solution which is 
robust against scalar disorder. Our predictions should be 
experimentally verifiable in HgTe QWs and Bi2Se3 thin 
films. 

The rest of this paper is organized as follows: In the next 
Section we give a short introduction to the model. In 
Section III we analyze the single QSHI/metal interface, 



which helps us to understand the metal/QSHI/metal 
junction results in Section [TV] We use tight-binding cal- 
culations in Section [V] to test the robustness of conduc- 
tance maximum in the presence of scalar disorder. We 
finish the paper with conclusions. 



II. BRIEF DESCRIPTION OF MODEL 

Close to the T point (k = 0), HgTe quantum wells or 
Bi2Se3 thin films can be described by an effective 4x4 
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Dirac-like modePC^nSI In the following we concentrate 
on the description of HgTe QWs while we will comment 
on Bi 2 Se 3 thin films at the end of this section. In the 
absence of structure inversion asymmetry the Hamilto- 
nian for HgTe QWs consists of two decoupled blocks, 
which are related by time reversal symmetry (TRS). Us- 
ing the basis (|-E1+) , \H1+) , \E1-) , \H1-)) T , with (El) 
and (HI) designating electron and heavy hole sub-bands 
and ± denoting the Kramers' partner, the Hamiltonian 

reads (k± 
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with d a (k\\) = (Ak x ,—Ak y ,AA(k)) T . a a are the Pauli 
matrices and I is the unit matrix in the subband 
space. Here e(fc) = C — Dk 2 is the energy dispersion, 
A4(k) — M — Bk 2 is the k dependent gap and A gives 
the strength of the coupling between electron and heavy 
hole states. The parameters A, B, D and M depend 
on the quantum well width. Below (above) the critical 
thickness of the quantum well M > [M < 0), the 
system describes the topologically trivial (non-trivial) 
regime. Due to the block diagonal form of Eq. (fTl) 
we will just consider the upper block in the following. 
The results for the lower block then follow by TRS. 
A generalization including a finite value of structure 
inversion asymmetry was presented in Ref. 1151 



In the following we consider the transport along x- 
direction through a single interface between QSHI and 
normal metal as well as QSHI connected to two metallic 
leads (see inset to Fig. [I] and Fig. 3a respectively) . We 
set C to zero in the center region, so called quantum spin- 
Hall insulator regime, and to a large negative value in the 
leads, thereby effectively modeling the metal contacts by 
highly doped HgTe, cf. Fig. [3a] for a schematic. So far 
similar systems were described within periodic boundary 
conditions (PBC^P^ where due to the periodicity in y- 



direction k y is conserved and quantized i.e. k y — 2-kti/W , 
where n — 0, ±1, . . . , ±i\r roax with iV max given by the 
number of propagating modes. However PBC do not 
allow for a formation of edge states and the competition 
between edge and bulk transport was not explored so far. 
Therefore in this paper, in order to include edge states 
we choose hard wall boundary conditions (HBC) i.e. the 
wave function is zero at the edges of the sample. An an- 
alytical solution of Eq. ([I]) for HBCP^ couples different 
k y modes and a simple wave matching procedure like e.g. 
in Refs. [TS] and [TH is no longer possible. To resolve this 
issue we expand the states in terms of the Fourier modes 
^n(y) = sin [mry/W], whose orthogonality allow 

an independent matching for each mode. Therefore, for 
an infinite quasi-lD system with k x conserved, the eigen- 



vector can be expanded as follows: 

oo 

i> m (x,y) = exp [ik?x] XnK (y) (2) 

n=l 

Here the index m labels the /c™-eigenvalues, Xn ^ s the 
two component spinor corresponding to the upper block 
of Hamiltonian (1) for a given transverse mode n. To find 
the eigenvalues and eigenspinors we use the Schrddinger 
equation 

oo 

(H - E) exp [ik^x] £ XZ <t>m (y) = 0. (3) 

Multiplying from the left by &y4>n 2 (y) an d using the 
orthogonality of the sine functions yields: 
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where all terms of the Hamiltonian proportional to k y — 
—id y are collected in H ky , all constant terms in % const 
and all with k x (k%) in U k * (U k *). Additionally we 
introduced x'™ — ^™X™- Defining the vectors x™ = 
(X™=1,X™ =2 ,---) T and X ' m = (x^X^.-.r, which 
are built by the 2 component spinors x™ an d X^T> Eq. Q 
can be written as matrix equation 
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So doubling the dimension of the system of equations by 
introducing x' m allows to reduce the p roblem of finding 
k™ to a linear eigenvalue equatio n 1 19 * 20 !. The sub-matrices 
in Eq. ^ 
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Further we used here D± = (D ± B) and E± = C - E ± 
M. The only term coupling different modes is rj nin2 = 
((j) ni | k y \(j)n 2 )- 

The solutions fc™ can be characterized as propagating 
(k x l e R) or evanescent (Im(fc x ) ^ 0). For real fc™ the 
propagation direction can be determined by the sign of 



dy^L(x,y) [dk x h(k)}u^ km Tp m {x,y). (10) 
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Analogously evanescent states with lm(k x ) > 
(Im(fc x ) < 0) are decaying to the right (left). In the 
following we will denote right (left) going solutions by 
the index run {mi,)- In general the Fourier ansatz also 
works when the two blocks of Hamiltonian ([!]) are con- 
nected by the Rashba-spin orbit coupling like it was e.g. 
done in Ref . QH 

Having the solutions for a single system we can now 
consider QSHI/metal and metal/QSHI/metal junctions. 
The interfaces in Figs . [T] and [3a| break translational sym- 
metry in x-direction and an incoming state ip- ln can be 
scattered in all possible k x -modes. We will now label 
the k x eigenvalues as: k x n 1 for the part 1 stretching from 
x = -co to x = -L/2, k™ 2 for -L/2 < x < L/2 and fc™ 3 
for x > L/2. If the incoming state is assumed as incident 
from the left in mode m, the state in the left lead takes 
the form 

ip mL ,i(x,y) (ii) 

and that of the right lead is 
^z{x, y) — E 

m R ,z{x,y). (12) 

rriR 

Here \r mL ^ m \ 2 and |i mji ,m| 2 are the probabilities that the 
incoming mode m is reflected into mode of the left 
lead or transmitted into mode m# of the right lead, re- 
spectively. Only in the central region of length L there 
are left and right outgoing states 

&2(x, V)=y~] c mRim ip mRt2 (x, y) + E d mL>m ip mLt 2(x, y). 

m R m L 

(13) 

The scattering amplitudes are computed by match- 
ing the wave functions ^i{x,y) and the associ- 
ated currents Jf{x,y) = [dk x h(k)] k _^_ i9 ^i(x,y) at 
the interfaces, i.e. * x (-L/2,y) " = § 2 (-i/2,y), 
9 2 {L/2 t y) = ML/2,y), J?(-L/2,y) = Jf(-L/2,y) 
and Jf(L/2,y) = Jg(L/2,y). For setups like in Fig. [TJ 
where we have only one interface (L = 0), the matching 
simplifies to *i(0,») = * 3 (0,y) and Jf(0,y) = Jf (0,y). 
The total transmission then is 

T= E E ^K R ,n R \ 2 . (14) 

Moreover the knowledge of all scattering amplitudes 
determines the full state ^(x,y) up to normaliza- 
tion. ty(x,y) is the response to an incoming mode 
i/jm,i(x,y) from the left lead. Therefore it allows 
us to calculate the non-equilibrium charge density 
n(x,y) = \^>{x,y)\ 2 as well as the non-equilibrium 
current density, J(x,y) — & (x,y)\7kh(k)^(x,y). 
For actual computation the Fourier series has to be 
truncated at a sufficiently high mode N. The error 
of this approximation depends weakly on the width 



of the sample and on the Fermi energy. However, 
the asymmetric shape of an edge state is harder to 
approximate by sine functions than bulk states. Nev- 
ertheless a good approximation can be achieved. For 
example the maximal relative errors in the 99 lowest 
fc™ eigenvalues at W — 1000 nm and Ef = changes 
from 2 % at TV = 100 to 0.6 % at TV = 200. The 
higher modes decay very fast and have little influence 
on the transport. Unless otherwise specified we will 
use the following parameters: A = 0.375 nmeV, B = 
-1.120 nm 2 eV, D = -0.730 nm 2 eV, M = -3.10 meV. 
When we are referring to the normal regime we put 
M = +3.10 meV. 

In Bi 2 Se3 thin films the overlap of surface states of 
oppo site s urfaces leads to the opening of a gap at V 
poiirpH22] This gap closes as the energetically lowest 
lying states, which have opposite parity, cross at cer- 
tain thicknesses d™ (n = I, 2, . . .) of the thin filnPI. At 
each crossing the parity changes leading to the oscilla- 
tions between topologically trivial and non-trivial insu- 
lators as a function of the film thickness. Around the 
band crossing points a foldin g of the Hamiltonian of the 
3D topological insulato r^ 23 * 24 ! to the l owest en ergetic sub- 
bands reproduces the model Eq. f or the 2D 
topological insulator. We choose for the calculations 
d z = 3 nm with the parameters of 2D effective model 
as follows: A = 0.406 nmeV, B = -0.250 nm 2 eV, D 
I). 070 nm 2 eV, M = -22.5 meV 1 ^. 



III. SINGLE INTERFACE BETWEEN QSHI 
AND METAL 

Here we analyze the injection of the helical edge states 
from the quantum spin Hall insulator into a metallic lead. 
The latter is modeled by Eq. ([T]) with a high doping, i.e. 
large |C|. The used setup is sketched in the inset of Fig. [T] 
On the QSHI side, the Fermi energy is chosen to be zero 
while on the metal side Ef is shifted by C, with C char- 
acterizing an energetic height (strength) of the interface. 
In Fig. [T] we show the transmission through such an inter- 
face as a function of C for different widths of the system 
W = 250, 1000 and 2000 nm. While for narrow sam- 
ples strong interfaces introduce a significant backscatter- 
ing due to an overlap of counter-propagating edge states, 
perfect transmission is observed for wide samples. The 
perfect transmission for wide samples can be explained 
analogously as for the topologically protected edge states 
in graphene^, where it was found that the counter prop- 
agating edge states for the same Kramers' partner are 
orthogonal. Indeed taking into account an exponential 
decay of the wave function for an edge stated it can 
be easily confirmed that the overlap between edge states 
goes to zero for W — > oo in the presence of any x depen- 
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FIG. 1: The conductance as a function of the energetic 
height of the interface at a junction between a QSHI 
and a doped QSHI. A sketch of the system is shown in 

the inset. Perfect transmission is achieved in wide 
samples where counter propagating edge states for the 
same Kramers' partner (spin) do not couple. Negative 
(positive) C indicate a QSHI/n-type metal 
(QSHI/p-type metal) interface. 



dent potential: 

(* k jV(x) - V(x) / dy*l*. kx W -^° 0. 

Jo 

(15) 

By the Onsager- Casimir symmetry relation^S122l injecting 
electrons from the metal to QSHI should be also perfect 
for very wide samples. Let us mention for completeness, 
that a linear Rashba coupling^ does not increase the 
backscattering, because, although the direction of edge 
state polarization is rotated, there are still two orthogo- 
nal solutions for edge states at the same edge. 



IV. METAL/QSHI /METAL JUNCTION 

In this section we consider metal/QSHI/metal junc- 
tions with a finite length L of the QSHI as shown in 
Fig. |3a| So far such junctions have been only analyzed 
using PBC which neglect the topological edge states^. 
In particular it was found that such junctions allow for 
the distinction of different topological phases purely due 
to evanescent bulk modes 12 . While the conductance rises 
monotonically with decreasing L in the topologically triv- 
ial regime (M > 0) , the topologically non-trivial regime 
(M < 0) is characterized by a conductance maximum at 
L max • In the limit of D = and BC > A 2 > BM 
(C < 0) the position of the maximum can be predicted 
for the zero mode by 

A , (2BC\ 
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FIG. 2: The conductance through a metal/QSHI/metal 
junction, like depicted in Fig. |3a[ as a function of the 
sample length L for different W (all values in nm). We 
compare the normal regime [M = 3.10 meV > 0) with 
the inverted one (M = -3.10 meV < 0) for HBC and 
PBC. Metallic leads are treated as highly doped HgTe 
with C = —2.5 eV. The inset shows the corresponding 
conductance signals for a Bi2Se3 thin film 
(M = ±22.5 meV) of width 250 nm at C = -2.5 eV. 



The purpose of this paper is to better understand the ex- 
istence of the conductance maximum in the topologically 
non-trivial regime as well as to study the interplay of the 
bulk and edge contributions using hard wall boundary 
conditions. We use the technique described in Ref . IT21 to 
solve the problem with PBC and the method described 
in Section |H] for a full solution with the hard wall bound- 
ary conditions. In Fig. [2] we present the conductance 
as a function of the sample length L. First we observe 
that the conductance decreases exponentially with the 
increase of L in the normal regime, independent of the 
boundary conditions. Indeed in this regime we expect 
only evanescent solutions and therefore the same results 
for PBC and HBC. We now focus on the topological in- 
sulator regime i.e. inverted regime with M < 0. In this 
regime, for HBC, both evanescent and edge states con- 
tributions to the conductance are present and the signal 
depends on the length and the width of the sample. 

In the very narrow QSHI regime, see black line in 
Fig. [2j we find Fabry-Perot like oscillations in the trans- 
mission of the edge states for HBC. At large L, the differ- 
ences between two oscillation maxima is given by 7r/fc c d gc , 
where fc c dgc is real wave vector of the edge states. The 
Fabry-Perot oscillations originates from the backscatter- 
ing of edge states at the interfaces and therefore are only 
present for narrow samples with the finite overlap be- 
tween counter propagating edge states. In this case the 
existence of edge states has a dramatic effect on the be- 
havior of conductance as a function of L. The quantum 
confinement in narrow wires shifts the energy levels of 
the states to higher values. This also increases the cor- 
responding imaginary k x values, which leads to a faster 
decay of incoming electrons in the QSHI tunneling junc- 
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tion. The combination of the Fabry Perot oscillations 
of the edge states and the bulk modes changes the posi- 
tion and the shape of the conductance maximum, which 
does not coincide with the maximum generated by cither 
the evanescent bulk states or the Fabry Perot oscillations 
alone. Therefore Eq. ( 16 1 can not be used to find a max- 
imum of conductance. 

In contrast, for a wide and long QSHI part, the conduc- 
tance is dominated by the edge states contribution and 
equals 2e 2 /h. For wide and medium length of QSHI, 
PBC and HBC give very different results while for short 
and wide QSHI the transport is dominated by tunnel- 
ing evanescent modea^. This is a very interesting result 
since naively one would expect that the edge state and 
bulk contributions to the conductance are additive for 
wide QSHI. 

The same behavior can be found for Bi2Se3 thin films. In 
the inset of Fig. [2] we show the conductance for a 250 nm 
wide strip. Instead of a maximum we observe a plateau 
like behavior in the inverted regime. However, increasing 
further the C parameter restores the conductance maxi- 
mum. Despite the narrow width the signal of HBC and 
PBC coincides around the plateau. This can be explained 
by the larger gap in Bi2Se3, which decreases the overlap 
of the edge states^. 

Since we showed in the last section that the transmission 
through a sufficiently wide metal-QSHI interface is per- 
fectly quantized, the scattering from the interfaces can- 
not be the reason for the absence of the edge state con- 
tribution. Therefore to understand the non-additive be- 
havior of bulk and edge contributions and the formation 
of the conductance maximum, we plot in Fig. [4] the evo- 
lution of the charge density for different sample lengths. 
In Fig. [4] we take into account only scattering of the first 
incoming mode with n=l, like it was described in Sec. |TT] 
after Eq. ( fl4| ), and use again the parameters for HgTe 
QWs. Since the lowest modes give the largest contribu- 
tion to the conductance maximum as shown in Fig. |3b[ it 
is justified to analyze first incoming mode signal to grasp 
the important details of the conductance maximum. For 
plotting the charge density we use the following normal- 
ization: 



tr 



dy& (x 
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Let us first concentrate on the Figs . [4a[ [4c| and [4d| where 
L 7^ L max . Due to the large C parameter a strong reflec- 
tion occurs at the left interface leading to standing waves 
in the left lead, see e.g. Fig. |4d| Further we observe that 
the density is peaked inside the sample. The incoming 
states from the left lead are purely metallic bulk ones 
since for high energies (high C parameter) there is no in- 
fluence of the edge states. In the short QSHI part (see 
Fig. 4a), these metallic states are not able to adjust to 



the shape of edge states and so they keep their form and 
tunnel through the sample. Therefore for short L the 
conductance signals from HBC and PBC should coincide 
as indeed found in Fig. [2] . When the length of QSHI 
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FIG. 3: (a) A schematic of the setup, consisting of the 
two metallic leads and the central QSHI at zero Fermi 
energy. The energy shift C is also shown in the band 
structures of the individual segments, (b) The 
conductance as a function of L for single modes 
(W = 1000 nm and C = —2.5 cV). The maximum is 
mainly carried by the low modes. 



part is comparable to the decay length of the evanescent 
modes, (see Fig. 4c), electrons are more likely transmit- 
ted through the sample by scattering from the bulk states 
into the edge states. For longer QSHI parts (see Fig. 4d) 
transport is driven by edge states. 
It is interesting to compare Figs 



4a 4c and 4d with 



Fig. |4b[ where we choose the length of QSHI part corre- 
sponding to the peak in conductance i.e. L « L max . One 
can see that the charge density is much larger in the QSHI 
part in comparison with other cases where L =/= L max . 
Additionally we do not observe standing waves in the 
left lead, which indicates a weaker interface reflection. It 
also becomes clear that at the position of the maximum 
the edge states are not yet populated, which is the rea- 
son for a coincidence of conductance signals for PBC and 
HBC. 

To further investigate the origin of the conductance max- 
imum, when bulk states dominate transport, we analyze 
the effective wave vector in the direction of propagation 
for the QSHI part 

W 

kf(x) = [ dy^Ux, y) (-icy * 2 (z, y), (18) 



where ^2{x,y) is the full scattering solution from 
Eq. ( |13[ ). Our analysis will focus on the lowest modes 
in PBC and HBC, which give the largest contribution to 
the signal. In Fig. [5] we present the imaginary as well as 
the real (inset) part of k° (x) for different lengths and 
the lowest incoming mode. The lowest modes for PBC 
and HBC are n=0 and n=l, respectively. In general fc° 
can exhibit a finite real part even for PBC, i.e. in ab- 
sence of edge states, due to evanescent solutions com- 
bined with complex amplitudes in ^2(x,y). Indeed one 
can see in Fig. that the real part of k% s is non-zero and 
decays exponentially away from the two interfaces. The 
imaginary part instead has a more sophisticated behav- 
ior with several crossings of zero. In particular Imfc° ff (x) 
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FIG. 4: Evolution of charge density in a 
met al/Q SHI/ metal junction as a function of the length 
of the sample. The charge response is the response to 
the first incoming mode, cf. Sec. Rafter Eq. (li]). The 
interfaces at ±L/2 are indicated by two black lines. In 

the plot, only upper edge state is seen since we 
restricted our calculations to the one spin direction. We 
used W = 1000 nm and C = -2.5 eV. (a) L = 50 nm, 
(b) L = 100 nm w L max , (c) L = 200 nm and (d) 
L = 500 nm. 



shows an antisymmetric behavior in respect to the sample 
middle for the length corresponding to the conductance 
maximum i.e. L = 100 nm w L max for lowest modes 
in PBC and HBC designated by blue/dark gray dashed 
and red/gray solid lines, respectively. This suggests that 
the imaginary part of vanishes after averaging over 
x-direction. We define (fc° ff ) as: 



L/2 

(kf)= J dxkf 

-L/2 



(x) 



(19) 



Fig. [6] shows (fc° ff ) as a function of the insulator length 
for PBC. The dashed red/dark gray curve in Fig.[6]shows 
that Im(fc° ff ) vanishes indeed at L max for the zero mode. 
The vanishing evanescent part of the effective wave vec- 
tor coincides with a maximum of the real part of (k° s ) 
(solid red/dark gray line), which gives the impression of 
an effectively propagating state leading to the conduc- 
tance maximum. Away from £ max , Ke(k° s ) decreases 
while Im(fcj S ) has a finite value for the zero mode. Indeed 
Fig. [5] shows that the signal for Imfc° ff and L = 2 * L max 
is not antisymmetric. The inset to Fig. [6] shows that the 
higher modes contribute less to the conductance maxi- 
mum, but for n = 1 there is still a Re(fc° ff ) maximum 




FIG. 5: The imaginary and real (inset) part of fc° ff 
value as a function of the position x for different lengths 
of QSHI part, boundary conditions and incoming modes 
n. We used W = 1000 nm and C = -2.5 eV. The 
conductance maximum corresponds to the length 
imax ~ 100 nm. 



around L max (solid green/light gray line in Fig. [6]). This 
maximum for n = 1 is smaller in comparison to n — and 
is shifted the same way as respective conductance max- 
ima for n = and n = 1 (see inset to Fig. |6j) . Im(k° s ) 
for n = 1 has also a zero value at £ max . In contrary for 
n = 5, which does not contribute to the maximum, the 
Im(fc° ff ) (black dashed line) is constant and the Re(k° s ) 
(black solid line) drops to zero. Therefore higher modes 
( modes with n > 4 for our set of parameters) behave al- 
most like modes in the topologically-trivial regime (com- 
pare rose/light gray thick lines with black lines in Fig. [6]) 
and have evanescent character. 

Another way to find differences between the topo- 
logically trivial and non-trivial regimes is by analyzing 
the conductance as a function of the Fermi energy for a 
sample of fixed length. For PBC the two regimes show 
the same behavior as long as the Fermi energy in the 
sample lies within conduction or valence bancP^l Only 
within the gap the QSHI can be distinguished from the 
normal insulator by a conductance peak around Ef = 0. 
For HBC we find the same behavior as within PBC in 
regimes where the bulk states dominate the transport, 
i.e. in wide and short QSHI junctions. When we increase 
the length of the QSHI above the decay length of the 
evanescent modes, the conductance is governed by the 
edge state contribution. In this limit within PBC the 
signal goes to zero while for HBC the conductance is 
2e 2 /h in wide wires where Fabry-Perot oscillations are 
negligible. In narrow and short samples however the 
transport is mediated by a mixture of edge and bulk 
states. Exactly this latter limit is shown in Fig. [7] 
where the Fermi energy dependence of the conductance 
is presented. In Fig. [7] , the bulk gap is indicated by 
two black vertical lines. It is larger than 2\M\ due to 
quantum confinement. Additionally the dashed vertical 
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FIG. 6: (fc£ ff ) as a function of the length for different 
incoming modes. Solid and dashed lines show real and 
imaginary parts of (k° s ), respectively. For clarity, only 
the calculations for PBC are shown. We see that for low 
modes the position of the conductance maximum (see 
inset) coincides with a maximum in the real part and a 
zero of the imaginary part of (k^ s ). For higher modes 

this behavior gradually disappears. We used 
W — 1000 nm and C = —2.5 eV. The maximum lies at 
imax ~ 100 nm. 
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FIG. 7: The conductance as a function of the Fermi 
energy in the central part of the junction. The Fermi 
energy in the leads is C = —2.5 eV. The sample is 
250 nm wide and 100 nm long. The vertical lines 
indicate the bulk gap (solid) and the mini-gap (dashed). 
The inset shows the density of states for spin-up 
electrons as a function of the Fermi energy and the 
y-coordinate in an infinite wire and for the inverted 
regime. The mini-gap can be identified by the lack of 
any propagating solution. 



lines present the energy range, in which the overlap of 
the edge states produces a gap in the system^. We will 
refer to it as mini-gap from here on. One can notice 
that outside of the bulk gap (solid lines) the signals of 
the normal and the inverted regi mes b ehave similar, like 
it was found for wide structure d 12 * 13 !. Within the gap 
however and for the PBC, a clear conductance peak in 



the inverted regime around Ef — is seen, while in 
the normal regime the transport is almost completely 
suppressed. For HBC one observes a similar behavior 
for the normal regime, but the existence of the edge 
states changes the signal drastically for M < 0. Within 
the gap and for energies smaller than zero, we find the 
quantized conductance of the edge states. Around zero 
and for positive energies the signal is strongly reduced 
and resembles that of the bulk modes. This behavior 
can be understood by looking at the density of states 
in an infinite wire in the inverted regime. In the inset 
to Fig. [7] we show the corresponding \ip\ 2 as a function 
of the y-coordinate and the Fermi energy. As before we 
consider just one of the blocks, so that only one edge 
state (spin-up) is shown. The mini-gap as well as the 
bulk gap can be easily seen in this contour plot. We 
find that the edge state spreads over the whole sample 
shortly before the mini-gap opens. We have seen that 
the overlap of the edge state crucially influences the 
transport in metal/QSHI/metal junctions. Therefore let 
us note that the penetration dep th of th e edge state is in 
general a function of the energy^* 30 1 3 ^ like it is shown 
in the inset of Fig. [7] For D finite and smaller than 
zero this behavior is asymmetric. For negative energies 
the edge state is strongly localized at the upper sample 
border. The small penetration depth leads to the fact 
that the edge state exists also at energies outside of the 
energy gap. For our parameters it merges to the valence 
band at energies of about Ef — —80 meV. Up to this 
energy the strong localization ensures a good quantiza- 
tion which we observe in the conductance. In contrast, 
on the other side of the mini-gap (Ep > 3 meV) the edge 
state is widely spread, so that a strong backscattering 
occurs. Moreover it merges quickly to the bulk states 
around Ef — 10 meV. Therefore we believe that at these 
energies the transport is dominated by evanescent bulk 
states. Hence in this regime the conductance resembles 
the shape of the PBC signal. 



V. DISORDER 

In this section we want to analyze the robustness 
of the conductance signal in the presence of scalar 
disorder. Disorder breaks the translational invariancc 
of the sample which makes the application of the 
wave matching routine impossible. We go around this 
problem by using tight-binding calculations within the 
Landauer-Biittiker formalis m 32 '' 3 ' In the tight-binding 
approach we discretize the sample using a finite grid 
of say N x x N y points while the hard wall boundary 
conditions are naturally implemented by missing bonds 
at the borders. In principle we could have performed 
the whole analysis, which we reported in this paper, 
using tight-binding calculations. However, to check 
the validity of the conclusions for the experimentally 
relevant setups one needs to simulate large structures. 
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Further for a good approximation within the tight- 
binding model one needs a small lattice constant a^, 
i.e. a small distance between neighboring sites, because 
many fast decaying evanescent modes have to be taken 
into account. As a consequence we would have to use 
large systems with up to (2000 x 200 x 2) 2 (width x 
length x degrees of freedom) lattice points, which need 
extraordinary amount of memory and CPU time. 

To analyze the influence of disorder let us here vary 
the sample structure and the parameters, similarly 
as it was performed for PBC in Ref. H31 First we 
decrease the width of the sample to 150 nm. For PBC 
this reduction should not change the position of the 
maximum. However, for HBC the overlap of the edge 
state wave functions influences the signal in narrow wires 
drastically. The overlap can be reduced by increasing 
the bulk gap, i.e. the M parameter as it was shown in 
Ref. [T7| Further Eq. ( 16 ) shows that a larger M param- 
eter leads to a smaller L max . Therefore, in the following 
by increasing M by the factor of 10 to M — —31.0 meV, 
we expect the maximum around 10 nm instead of 
100 nm. Indeed Fig. [8] shows the conductance maximum 
at about 12 nm. First we compare the tight-binding 
(TB) calculations in the absence of disorder (red/gray 
dashed line in Fig. [8]) with the signal computed in the 
wave matching (WM) formalism (red/gray solid line). 
The signals of the two models almost coincide when we 
choose N y — 250 grid points. The number of points 
in x-direction varies with L, since we fix the lattice 
constant to = 150 nm/251. Therefore indeed the 
above described procedure does not change the results 
quantitatively but shifts the position of the conductance 
maximum to system sizes which can be easily calculated 
within tight binding model. Further the conductance for 
PBC (red/gray dotted line) coincide with the signal from 
tight-binding calculations indicating that the overlap 
between the edge states is small. 

Scalar disorder is included by choosing randomly 
Nd = UdN x N y (rid G [0, 1]) lattice sites. For our analysis 
here we fixed = 0.1. On each of the selected points 
we shift the energy randomly by Vd € [— W<j/2, Wd/2]. 
The strongest disorder we can put without coupling the 
valence to the conduction band therefore is Wd = 2\M\. 
The disorder is modeled by a Gaussian-like function 
with the decay length given in units of the lattice 
constant a^. All signals including finite disorder have 
been averaged over 100 different disorder realization. 
In Fig. [8] we focus on the constant disorder strength 
Wd = 0.06 eV varying the decaying length of the 
Gaussian disorder from rd = la^ (solid black line) to 
2ol (dotted black line). For long sample lengths L, the 
disorder has no influence on the quantized conductance 
of the edge states, because the scalar disorder cannot 
couple the counter propagating edge channels. When 
evanescent waves contribute to the transport at smaller 
L the disorder has a visible influence on the signal. 
The influence of the disorder increases with an increase 




WM HBC ■ 
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TB, n d = 0. r d = . 

TB, n d = 0.1, r d - 1.0 a L ■ 

TB, n d = 0.1, r d = 1.5 a L ■ 
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FIG. 8: The charge conductance as a function of the 
sample length. We compare the wave matching (WM) 
approach for HBC (solid red/gray line) and for PBC 
(dotted red/gray line) with the tight-binding 
calculations (dashed red/gray line) and find good 
coincidence. Adding scalar disorder with Gaussian 
correlations influences the signal slightly but does not 
destroy the conductance peak. We used W = 150 nm, 
C ~ —2.5 eV and M = 31 meV. The disorder strength 
was set to Wd — 0.06 eV for the black lines. The inset 
shows the influence of the Fermi energy Ef inside the 
sample region. 



of disorder decay length. It influences the signal by 
shifting and increasing the maximum of conductance. 
To understand this effect we need to remember that 
the transport for short QSHI lengths is dominated 
by evanescent modes. In this regime the extension of 
impurities can be comparable with the sample length 
introducing large fluctuations of energies over a relatively 
large part of the sample which can cause a shift of the 
conductance maximum . Indeed the inset to Fig.|8]shows 
the conductance signal shift to left or right depending on 
the Fermi energy inside the QSHI. Despite these small 
changes in the conductance behavior, the conductance 
signal is robust against disorder for HBC. 



VI. CONCLUSION 

We analyzed the interplay of edge and bulk states in 
metal/quantum spin-Hall insulator (QSHI)/mctal junc- 
tions as a function of the size of QSHI part as well as the 
Fermi energy. We found that for short and wide QSHI 
part, the edge states are not populated due to conduc- 
tance mismatch between highly doped leads and QSHI 
regime, and the transport occurs through the evanescent 
modes. In this regime, the topologically non-trivial in- 
sulator can be distinguished from the trivial insulator 
by the appearance of the distinct conductance maximum 
at small QSHI length and the conductance peak in the 
gap when the Fermi energy changes for the constant 
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length. Moreover, we showed that the origin of this con- 
ductance maximum is surprisingly the formation of an 
effectively propagating state with the real k x vector built 
from evanescent modes. 

In contrast, in the wide and long samples the transport 
is driven by spin edge channels while in the intermediate 
size regime the conductance signal is built from both bulk 
and edge states in the ratio which depends on the width 
of the sample. In this latter regime the combination of 
the Fabry-Perot oscillations (coming from the overlap of 
edge states) and bulk modes changes the position and 
the shape of the maximum. 

Our predictions are robust against scalar disorder. 
Therefore the measurements of systems with diferent geo- 



metrical aspect ratios should allow to distinguish the bulk 
and edge state contributions to the transport. Further 
the detection of the bulk conductance peak as a function 
of Fermi energy or the length for short and wide QSHI 
junctions, should be an alternative method to distinguish 
between topologically trivial and non-trivial insulators in 
HgTe QWs as well as Bi 2 Se3 thin films. 
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